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Abstract 

A billiard in the form of a stadium with periodically perturbed boundary 
is considered. Two types of such billiards are studied: stadium with strong 
chaotic properties and a near-rectangle billiard. Phase portraits of such bil- 
liards are investigated. In the phase plane areas corresponding to decrease 
and increase of the velocity of billiard particles are found. Average velocities 
of the particle ensemble as functions of the number of collisions are obtained. 
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1 Introduction 

A notion of billiard in physics is known since G.Bikhoff [0 who considered a prob- 
lem concerning the free motion of a point particle (billiard ball) in some bounded 
manifold. So, the billiard dynamical system can be introduced as follows. A billiard 
table Q is a Riemannian manifold M with a piecewise smooth boundary dQ. The 
billiard particle freely moves in Q. Reaching the boundary, it is reflected from it 
elastically. Thus, the billiard particle moves along geodesic lines with a constant 
velocity. In the present article we consider billiards in Euclidean plane. In this case 
the angle of incidence of the particle is always equal to the angle of reflection. 

In accordance with the boundary geometry, dynamics of the billiard particle can 
be integrable [^], completely chaotic ^ and, depending on the initial conditions, 
regular or chaotic 0, |^, |, [l^ . 

A natural physical generalization of the billiard problem is perturbation of the 
boundary in one or another manner. For the first time, this problem concerning 
collisions of particles with massive moving scatterers has been considered by S.Ulam 



1^ in the context of the unbounded increase of energy in periodically forced Ha- 
miltonian systems. It goes back to the question related to the origin of high energy 
cosmic particles ||ll[ and known as Fermi acceleration. The Fermi-Ulam model has 



been the first system where invariant curves, chaotic layers and stable islands have 
been investigated (details see in [l^). It has been shown that in the case of a quite 
smooth perturbation of the boundary the particle velocity is bounded by invariant 
curves. Otherwise, the velocity can grow indefinitely. 

For low- dimensional integrable billiards the problem of Fermi acceleration has 
been studied on the example of circle and elliptic billiards |^, |14| . In these papers 
the authors come to conclusion that the velocity of the particle ensemble is bounded 
by the corresponding invariant curves. Investigations of chaotic billiards have been 



performed for the Lorentz gas O, |16|. As predicted, perturbations of the boundary 



of such a billiard lead to the appearance the Fermi acceleration for the particle. In 
addition, the acceleration is higher in the case of periodical boundary oscillations 
than in for their stochastic perturbations. 

In the present paper we study so-called stadium-like billiards [|lOl which are 
defined as a closed domain Q with the boundary dQ consisting of two parallel lines 
and two focusing curves (Fig.|l]). If parameter 6 is a sufficiently small then the 
billiard is a near-integrable system. In this case its fixed points are stable. As a 
result, in the stochastic (or chaotic) "sea" the stability regions appear which consist 
of invariant curves. At the same time, owing to a weak nonlinearity, dynamics 
near the separatrixes divided the stability regions of elliptic points is chaotic, and 
the particle can reach neighbourhoods of all points in the chaotic layer. In the 
case of the fixed boundary the particle dynamics can be both chaotic and regular, 
depending on the initial conditions. Introduction of external perturbations leads to 
the possibility of the particle passage from chaotic region to the regular one and 
back. This is the reason of new interesting effects which are also described in the 
paper. 



2 Definitions and maps 

In this section, basic analytical results are presented. They are necessary for the 
further description of the billiard dynamics. 

2.1 Stadium-like billiard with the fixed parabolic focusing 
components 

Consider a billiard shown in Fig.Q. To describe its dynamics let us construct the 
corresponding dynamical system for h <^ a. For this, we use the known method of 
specular reflections. It consists of the reflection of the billiard table form neutral 



components. As a result, the stadium is replaced by a caterpillar billiard. It can be 
shown that the change in the particle velocity in both cases is the same. In addi- 
tion, one can show that between trajectories of the initial billiard and the obtained 
"caterpillar", a one-to-one correspondence takes place. 




x„ a 

Figure 1: A stadium-like billiard and its development. 

Suppose that at the initial time the particle belongs to the billiard boundary and 
its velocity vector directs towards the interior of the billiard region. Let us choose 
coordinates ip and x as shown in Fig.|l|. The motion of the billiard particle induces 
a map (■?/'„, x„) (ipn+iyXn+i)- Suppose that h <^ I. In this case the focusing 
components can be approximated by the function = 4bx{x — a) /a?. For such 
billiard configuration the map is written as follows: 

Xn+i = Xn + l tan ipn+i , mod a , 

iJn+l =i'n- 2/?(x„+i) , 

where (3{x) = arctan (see Fig.|I|). If b is small then (3 ~ 4b{2x — a) /a?. For 

the further analysis change the variables: ^ = x/a, G [0, 1). Then 

/ 

in+i = in + - tan %l)n , mod 1 , 
a 

(1) 

86 

■0,1+1 = (2^n+l - 1) • 

a 

In Fig.^the phase portrait generated by the map (|1|) is shown. Initial conditions 
for each trajectory are marked by crosses. One can see that the trajectory with 
initial conditions in the chaotic region can reach any point of this region. At the 
same time, in the regular regions the points moves along invariant curves. 
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Figure 2: Phase portrait of the biUiard with parabohc focusing components (see map (0j ) 
at a = 0.5, b = 0.01 and I = 1. The diagram contains three regular trajectories (each by 
10^ iterations) and one chaotic trajectory (5 ■ 10® iterations). 



It is obvious that the fixed points of the map are the following: ^ = 1/2 and 
ips = arctan(ma//). If m = then the billiard particle moves strictly vertically. If 
m = 1 then it moves for one cell. And so on. In Fig.^ the fixed points correspond to 
maximal ellipses. Let us find the stability conditions of these points. To this end, 
change the variables: = + 1/2, ipn = ^'ipn + arctan(ma//) and linearize the 
map. Then we get: 



n+l 



a cos^ ip. 
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where ipg = arctan(ma//). The corresponding transformation matrix has the form: 

/ 1 -J— \ 



A 
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1 - 



a cos^ 
166/ 



cos^ ips I 



It is not hard to see that detA = 1. Thus, the map preserves the phase volume. 

The stabihty criterion for the fixed points is |TrA| < 2. Then cos^-?/'s > Abl/a^ 
or < //(46) — P/a^. On the other hand, transition to chaos take place if 

^?>l. (2) 

Eigenvalues of the matrix A are Ai,2 = e^*'^, where cos a — ^Tr^l. Let us introduce 
f — //(acos^Vs)) 9 — i^b/a. In this case 

\-g 1-fg 

and its eigenvalues 

X,,2 = ^ e^'" - 1 

Consider matrix X with the columns of eigenvectors. As known, in this case the 
matrix A = X~^AX is diagonalization of the matrix A: 



A 



e" 



New variables for which the transformation matrix has a diagonal form are the 
following: 



Z* J \ Alp 

where 

Y-i = ' ( - 1 -f 

2sin(7 V + 1 / 

One can see that Z and Z* are complex conjugate. Thus, 

If we take Z in the form of Z = 7e*^ then in the action-angle variables we obtain: 

Thus, the particle motion around the stable point is written by the map with the 
following rotation number: 



In turn, old variables have the form: 



= 2Jcos^ , 
2/ 

Atp = — (cos(o- + 9) - cos 9) . 

4/ 

Therefore, the Jacobian of the transformation is J = — — sin a. 

2.2 Perturbations of the boundary and resonance 

The particle motion with the velocity V generates a flow for which we can introduce 
time t. In turn, the time between two sequential collisions is the following: r ~ 

. Thus, we get the rotation period: 

cos yjg V 

Ti = —T= ^ . 

^ cos'j/'s arccosi 1 — 86//(a cosV")^ I 



If the system undergoes external perturbations of the boundary of period Text 
then for Ti = T^xt we can observe the resonance between degrees of freedom. This 
leads to the fact that regions including stability areas (see Fig.|^) are accessible for 
the particle. 

Some words about the nature of the resonance (see Fig.^. At the motion along 
the invariant curve in the neighbourhood of the stable point, angle oscillate. 
Collision of the particle with the perturbed boundary leads to the change m. ip. If 
the boundary moves towards the particle, then the angle decreases. Otherwise it 
increases. Suppose that the image of the trajectory moves along the arc AB. In 
this case, if the particle undergoes collisions coming from the opposite side, then the 
trajectory tends to inside of the area. If the boundary and the particle velocities 
have the same direction, then the trajectory goes outside of the region. 

For the resonance between the boundary oscillations and the motion along the 
invariant curve form B to A, collisions change the particle direction. This leads to 
larger shift of the trajectory to the stable point. 

From the equality Ti = Text we can obtain the resonance condition for the particle 
velocity: 

Vr = 1^ V • (4) 

COS ips arccos I 1 — 8bl/{a cosips)"^ 
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Figure 3: Invariant curves around stable points (schematically). 

Hereafter, for simplicity we assume that the frequency of external perturbations 
LJ = 1. Then Text = 27r. 

For the invariant curves with ips = arctan(ma//), m > 1, there exists only one 
area where the absolute value of the angle ipn increases (along the arc AB from A to 
B in Fig.^), and the area of the decrease in ipn- However, for the central stable fixed 
point there are two such areas: CD, EF for increasing and DE, EC for decreasing. 
This leads to the fact that in the neighbourhood of this point resonance is observed 
for the same condition that in the other areas, but the particle velocity is less by 
half. Therefore, 

v° = — = ^ r^i 

2 arccos(l-86//a2) " ^ ^ 

2.3 Focusing components in the form of the circle arcs 

In this section, we consider a stadium-like billiard with the boundary consisting of 
two focusing components in the form of the circle arcs and compare the obtained 
results with the previous parabolic case. 

2.3.1 Fixed boundary 

Suppose that focusing components are arcs of the radius R circle (symmetrical about 
the vertical billiard axis) and the angle measure 2$ (Fig.^). Geometrically we can 
obtain that 

R = ; <P = arcsm — . 

86 ' 2R 

For such a billiard the chaoticity condition is obtained as follows 0. Assume that 

Q G R^ and the focusing component is a part of the circle C. Chaos can be observed 



if disk dD = C, belongs to the billiard table Q. Thus, 

J_ _ 46/ 

that is the same as obtained from the analysis of the stable points. 
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Figure 4: A stadium-like billiard with focusing components in the form of circle arcs. 

Let us introduce dynamical variables as shown in Fig.^. Assume that angles ifn 
and a* are counted counterclockwise, and the angle a„ is counted clockwise. For the 
fixed boundary a* = a;„. Suppose that Vn is the particle velocity, and t„ is a time of 
n-th collision. Let us find a map describing dynamics of the particle in such a billiard. 
For this it is necessary to consider two cases: 1) After collision with the focusing 
component the particle collides with the same component (multiple collisions); 2) 
After the collision, the particle moves to the opposite focusing component. 
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Figure 5: Multiple collisions with the focusing component. 



1) Multiple collisions (Fig.||). 



In this case, geometrically we get the following map: 



Otn+1 — , 

(fn+1 = fn + T^ - 2a„ (mod 27r) 
2_Rcos a„ 



(6) 



t 



n+l 



If Iv^n+il < then the particle collides with the same component. Otherwise, 
n + 1-th collision with the opposite components takes place. 

2) Collision with opposite components. 

For this case the map can be written as follows: 



a. 



n+l 



arcsm 



sin {ip„ + ^) 



X. 



n+l 

R 



COsV'ri 



V^n+l =1pn- "n+l , 

i?(cOS (fn + COS iPn+1 — 2 COS $) + / 



(7) 



t 



n+l 



Vn COS Ipn 



where 



Oir 



R 



[sin an + sin($ - 



Xn + ^ tan ipn (mod a) 



Really, let us extend the circle arc up to the semicircle (Fig.^. Introduce the 
angle ip between the vertical and the velocity vector. It is counted clockwise. It is 
obvious that BC = Rcos(pn+ita.nipn+i- From triangles ABC and ABiCi we get: 

BiCi ABi i?(cOS(y9„+i - cos$) 



BC 



AB 



-. Therefore, 



R cos Lfn+l 

BiCi = -Rtan?/'„+i(cos(^„+i — cos$) . 

Now, using the value Xn for the next collision we obtain: 

Xn+i = Rsin^ + R sin v^n+i + R tan ipn+i (cos v^n+i — cos $) 

R 



(8) 



COS^n 



(sina„ + sin($ - 



In addition, 

^n+i = Xn + Itaxitpn (mod a) . 

If we invert the particle motion then the expression (|^) gives connection between x* 
and a*. This fact allows us to get the explicit system (|^. 

2.3.2 Perturbed boundary 

Suppose that focusing components are perturbed in such a way that their velocity in 
each point is the same and directed by the normal to the component. Assume that 
the velocity value depends on time as follows: U{t) = Uofiuit + to)); where is a 
frequency oscillation. We consider the case Uojuj <^ I, i.e. the shift of the boundary 
is small enough and can be neglected. Therefore, the billiard map is written as 
follows: 

Vn 

an 



v/K^„i + 4K-icosa;;f/„ + 4^2 



.'^-1 . A (9) 
arcsm | — -— sm ' 

Vn 



a. 



Or 



V^n + TT — 2an (mod 27r) 
2R cos an 

yn 



(10) 



x'.' 



R 



[sin an + sin ($ - ?/'„)] 



n+l 

ctn+1 = arcsin 



Xn + Itamjjn (mod a) 



sin {ipn + $) - 



n+l 

R 



cos 



ifn+l ='ipn- Oin+l , 

i?(cOS + COS — 2 COS $) + Z 

tn+l 



tn + 



Vn COS Ipn 



if + TT - 2an\ > $ 



:ii^ 



The given map describes a stadium-like billiard with the focusing components in 
the form of the circle arcs. This map is exact one except for the approximation 
Uq/uj -C /. The first group ( p!OD corresponds to sequential multiple collisions with 
one of the focusing components, and the second group (0) corresponds to the 
passage to the opposite side of the boundary. 



3 Numerical analysis 



In this Section we consider stadium-like billiards with constant and perturbed bo- 
undaries. In the first case, the particle dynamics is described by the exact map (|l]) 
and approximate map ©-(l^), respectively. 

3.1 Phase diagrams of billiards with the fixed boundaries 

Phase portrait can help to understand the system dynamics and define chaotic and 
regular regions in its phase space. In Fig.|^ the phase diagram of the billiard with 
the fixed parabolic boundary is shown (see map (|l|). Crosses in this figure are initial 
conditions. One can see that phase plane is divided into regular and chaotic regions. 
If the initial conditions belong to the regular region then the trajectory remains 
here and forms the corresponding invariant curves. However, beginning in chaotic 
region the phase trajectory uniformly covers this region. The described portrait has 
been obtained on the basis of three regular trajectories (each contains 10^ iterations) 
and one chaotic trajectory (5 ■ 10* iterations) of the map (|I|). Geometric size of the 
billiard is the following: a = 0.5, b = 0.01, / = 1. 

In Fig.|^ similar results concerning the map (|§)-(0), where the depth of the 
focusing components is taken into account, are shown. Remind that for this map 
the approximation 6 -C Z, a has been used. To reasonable comparison, the geometric 
sizes of the billiard and the number of trajectories have been chosen the same as in 
the previous case (Fig.|^). 

The difference between obtained diagrams can be easily explain by means of 
Fig.0. 

Approximation of a small enough depth of the focusing component for the map 
(|l]) means the following. Let B be an intersection point of the parabola ends and 
the particle trajectory. In approximation, we consider the particle collision in the 
point A which is a projection of B into the arc. But, in fact, the point C is the 
collision point. As a result, in the exact case for large enough ip, collisions happen 
mainly with the right (in the Figure) part of the arc, and from collision to collision 
the angle ip decreases. Thus, the billiard particles as if "push out" to the region 
of the small angel ip. This is in agreement with Fig.^ where the region ip < 11/2 is 
empty.. 



Figure 6: Phase portrait of stadium-like billiard with focusing components in the form 
of circle arcs (see map Parameters of the billiard are the same as in FigJ|. One 

can see nonuniformity of the covering of the chaotic region. 
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3.2 Perturbed map 

In this section we consider the problem of the velocity change depending on its 
relation to the resonance value (^. 

3.2.1 Phase diagrams 

Consider the map (p|)-(pn]) corresponding to the perturbed stadium-like billiards. 
Construction of the phase diagrams have been performed for the same values of 
geometric parameters as in the previous §3.1. Therewith, the amplitude of oscil- 
lations Uo = 0.01. As noted above (see §2.2), for various particle velocities the 
corresponding phase portraits should be different form each other. 

In Fig.^the resonance velocity as a function of the angle ips (see (^)) is shown. 
One can see that in the region form to ips max (where ips max is a maximal angle for 
which the fixed points are still stable) the value of the resonance velocity is varied 
through a small range. 

Phase portraits for the perturbed billiard are shown in Fig.^. For detailed nume- 
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Figure 7: The difference between the exact hilliard map and approximate map. 
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Figure 8: The resonance velocity as a function of ips (see ^). 



rical analysis three particle ensembles with initial value Vq = 1, 1.2 and 1.5 have been 
considered. Therewith initial conditions have been chosen in the chaotic region in a 
random way. Thus, the obtained portraits give an insight into the velocity change 
of the billiard particle. In the obtained Figure, the vertical shaded areas correspond 
to the velocity increasing, and the horizontal shaded areas fit its decreasing. The 
wait areas (without shading) are the intermediate ones; here the particle velocity 
is transient. The black tones correspond to the areas which are inaccessible for the 
phase trajectory. 

As follows from obtained diagrams, if Vq 7^ K then around the stable fixed 
points there exist the areas surrounded by invariant curves. As before, these areas 
are inaccessible for particles from the chaotic regions. At the same time, in the 
neighbourhood which has become accessible for the particles as a result of pertur- 
bations, one can see areas of the increasing and decreasing velocity. Depending on 
the relation to the resonance velocity value, they can change places. 

U Vq = Vr (the resonance harmonic) then all neighbourhoods of the stable fixed 




Figure 9: Phase diagrams of the velocity change in the biUiard with the perturbed bound- 
ary (see (iJ-dlDJ at b = 0.01, a = 0.5, I = 1, Uq = 0.01 and lj = I. Vq = 1.0, 1.2 
(resonance), 1.5. 
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Figure 10: The same as in Fig.|| but Vq 



= 0.5, 0.6 and 0.7. 



points (except for the central one, ipo = 0, = 1/2) become accessible for the 
trajectory. In addition, for this resonance there are no the well-defined areas where 
the particles have an acceleration. 

Following d^) the resonance velocity in the neighbourhood of the central stable 
fixed points ~ 0.6 (primary resonance). In Fig.^ the phase diagrams for the 
initial velocities Vq = I5 1-2 and 1.5 are shown. One can see that at Vq = all 
areas of the phase space are accessible. 



3.2.2 The particle velocity as a function of iterations 



Numerical investigations of the perturbed billiard described by the map (|9|)-(p!l 
have been performed in two cases: when the billiard has strong chaotic properties 
and for a near-rectangle stadium. In the first case the billiard is a "classical" sta- 
dium. Then \E' = 7r/2 and the billiard is a domain with a boundary that consists of 
two semicircles and two parallel segments tangent to them. The latter case means 
that focusing components are segments of the almost straight line, and the billiard 
system is a near-integrable one. 

For the first case the following billiard parameters were chosen: a = 0.5, b = 0.25, 
/ = 1, Mo = 0.01, u = 1, and Vq = 0.1. The particle velocity was calculated as the 
averaged value of the ensemble of 5000 trajectories with different initial conditions 
(solid curve 1 in Fig.^). These initial conditions were different from each other 
by a random choice of the direction of the velocity vector vo- As follows from 
the numerical analysis, the obtained dependence has approximately the square-root 
behaviour {V{n) ~ y/n)- The fitting function y ~ an'^ (the dot-and-dash curve 1 in 
Fig.|Tl|) yields the following values: a = 0.01015 ± 0.00002 and c = 0.4446 ± 0.0002. 

A near-integrable case means that parameter b (see Fig.|^) is a sufficiently small, 
and the curvature of the focusing components gives rise only weak nonlinearity in 
the system. In such a configuration the billiard phase space has regions with regular 
and chaotic dynamics. This case is much more interesting for investigations. 

As follows from numerical investigations, on each side of the resonance the be- 
haviour of the particle velocity is essentially different. If the initial value Vq < Vc 
then the particle velocity decreases up to a certain quantity V/m < and the par- 
ticle distribution tends to the stationary one in the interval (0, V/j„). If, however, 
Vq > Vc then billiard particles can reach high velocities. In this case the particle dis- 
tribution is not stationary, and it grows infinitely. In addition, the average particle 
velocity is also not bounded. 

For detailed numerical investigations initial conditions were randomly chosen in 



the chaotic region of the unperturbed bilhard. In Fig.[TT| the particle velocity as a 
function of the number of iterations is shown (curves 2-5). The billiard parameters 
remains the same as for the stadium-like billiard (curve 1) except for b = 0.01. 
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Figure 11: Average velocity of the ensemble of 5000 particles in a stadium as a function 
of the number of collisions, I = 1, a = 0.5, Uq = 0.01 and uj = 1. Two (dot-dash and solid) 
curves 1 corresponds to the billiard with strong chaotic properties (b = 0.25). Curves 2-5 
correspond to the near-integrable system (b = O.Olj; Vq = 1 (curve 2, 3) and Vq = 2 (curve 
4, 5). Curves 2 and 4 are the average velocities of the particle ensemble. Curves 3 and 5 
correspond to maximal velocities reached by the particle ensemble to the n-th iteration. 

On the basis of 5000 realisations and for every initial velocity, three curves have 
been constructed: the average, minimal and maximal velocities which the particle 
ensemble has reached to the n-th iteration. So, we can see the interval of the velocity 
change. As follows from this Figure, if Vq < Vc then the averaged particle velocity 
(solid curve 2) gradually decreases and tends to a constant. The maximal velocity 
of particles (dotted curve 3) also decreases up to V/j„ and then fluctuate near this 
value. Eventually, the particle velocities lie in the interval < V < V/j„. In the case 
of > Vc, the minimal velocity of particles as before decreases. This means that in 
the ensemble there is a number of particles which are in the region of low velocity 
values. In our numerical analysis the part of such particles was about 75 percents. 
At the same time, there are particles with high velocities (dashed curve 5 which 
corresponds to the maximal velocity of the ensemble). As a result, the averaged 



particle velocity (solid curve 4) increases. 
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Figure 12: Stationary distribution of the particle velocity calculated hy the one particle 
trajectory during 10^ iterations. Vfm is a maximally reached velocity. 

In Fig.0 a stationary velocity distribution is shown. This distribution was cal- 
culated by the one particle trajectory during 10^ iterations. The initial velocity 
was chosen as follows: Vq ~ V/j„/2. The value denoted by Vfin corresponds to the 
maximally reached velocity. 

4 Concluding remarks 

Billiards are very convenient models of several physical systems. For example, parti- 
cle trajectories in billiards of specific configuration can be used in modelling a lot of 
dynamical systems. Moreover, most approaches to the problems of mixing in many- 
body systems go back to billiard-like questions. A natural physical generalisation of 
a billiard system is a billiard whose boundary is not fixed, but varies by a certain 
law. This is a quite new field which opens new prospects in studies of problems that 
have been known for a long time. 

In the present article we have studied the problem of the billiard ball dynamics in 
a stadium with the periodically perturbed boundary. Numerical analysis shows that 
for the case of the developed chaos, the dependence of the particle velocity on the 




number of collisions has the root character. At the same time, for a near-rectangle 
stadium an interesting phenomena is observed. Depending on the initial values, the 
particle ensemble can be accelerated, or its velocity can decrease up to the quite low 
magnitude. However, if the initial values do not belong to a chaotic layer then for 
quite high velocities the particle acceleration is not observed. 

Analytical description of the considered phenomena (mechanisms of deceleration 
and acceleration of the billiard particle, stabilisation of unstable points etc.) required 



more detailed analysis and will be published soon [1/ 
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